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Abstract 

We report improved measurements of B to pseudoscalar-vector decays containing an to meson 
in the final state. These results are obtained from a data sample that contains 386 million BB 
pairs collected at the T(45) resonance, with the Belle detector at the KEKB asymmetric energy 
e + e~ collider. We measure the branching fractions B(B ± -> u^) = (8.1 ± 0.6 ± 0.5) x 10~ 6 , 
B{B ± -» um*) = (7.0 ± 0.6 ± 0.5) x 10" 6 , and B(B° -» ujK°) = (3.9 ± 0.7 ± 0.4) x 10" 6 . We also 
set the 90% confidence level upper limit B(B° — > lot: ) < 1.5 x 10 -6 . In addition, we obtain the 
partial rate asymmetries A C p = 0.05 ± 0.08 ± 0.01 for B ± -► uK* and A C p = -0.03 ± 0.09 ± 0.02 
for B^ — > u>tt^~ . 



PACS numbers: PACS numbers: 13.20.H 



Charmless hadronic B decays provide a 
rich ground to understand the dynamics of 
B meson decays and the origin of CP vi- 
olation. Two-body B decays with a vec- 
tor meson u and a pseudoscalar particle h 



(h can be kaon or pion) proceed through 
b — > s penguin diagrams and b — > u diagrams. 
Theoretical expectations from QCD factor- 
ization approaches suggest that the branch- 
ing fractions of charged B decays are be- 
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tween 3 and 8 x 10~ 6 and that the uK + 
branching fraction is smaller than the ujtt + 
0, 0- The decay B° — > un is color sup- 
pressed while the neutral B decay to uK° is 
expected to have slightly smaller branching 
fraction than the charged B decays. Experi- 
mentally, clear signals have been observed in 
B ± -> ljK ± ,B ± -> uo^ and B° -> uK° 
with similar branching fractions 0, 0, II- 
However, experimental measurements are not 
yet precise enough to tell which mode has the 
larger branching fraction. Furthermore, the 
branching fraction of B° — > uK° measured 
by the Belle collaboration is smaller than the 
BaBar result. Therefore, it is necessary to 
update the results with more data. 

In this paper, we report improved mea- 
surements of branching fractions and par- 
tial rate asymmetries for B — > uoh decays, 
where h can be a kaon or pion. The partial 
rate asymmetry (Acp) is measured for the 
charged B decays and defined to be 

_ r(g-->/-)-r(i?+->/+) 

CP ~ T(B~ ^ f-)+T(B+ ^ f+y { ) 

where T(B~ — ► /~) is the B~ decay rate 
and T(B + — > / + ) denotes that of the charge 
conjugate mode. These measurements are 
based on a data sample of 386 million BB 
pairs, collected with the Belle detector at the 
KEKB asymmetric-energy e + e~ collider. 
The statistics are a factor of 5 larger than 
our previous published results. 

The Belle detector is a large-solid-angle 
magnetic spectrometer that consists of a sil- 
icon vertex detector (SVD), a 50-layer cen- 
tral drift chamber (CDC), an array of aero- 
gel threshold Cerenkov counters (ACC), a 
barrel-like arrangement of time-of-flight scin- 
tillation counters (TOF), and an electromag- 
netic calorimeter comprised of CsI(Tl) crys- 
tals (ECL) located inside a super-conducting 
solenoid coil that provides a 1.5 T magnetic 
field. An iron flux-return located outside of 
the coil is instrumented to detect K\ mesons 
and to identify muons (KLM). The detector 
is described in detail elsewhere JzJ- In Au- 
gust 2003, the three-layer SVD was replaced 



by a four-layer radiation tolerant device. The 
data sample for this analysis consists of 140 
fb of data with the old SVD (Set I) and 
217 fb -1 with the new one (Set II). 

Hadronic events are selected using crite- 
ria based on the charged track multiplicity 
and total visible energy with an efficiency 
greater than 99% for generic BB events. All 
primary charged tracks must satisfy quality 
requirements based on their impact param- 
eters relative to the run-dependent interac- 
tion point (IP). For tracks from the candidate 
B mesons, their deviations from the IP posi- 
tion are required to be within ±0.1 cm in the 
transverse direction and ±3. 0cm in the lon- 
gitudinal direction. Charged particle iden- 
tification is performed using the K-tt like- 
lihood ratio, TZk = C>k/{£>k + where 
£>k {£>■*) is the likelihood of a charged par- 
ticle to be a kaon (pion) based on the infor- 
mation coming from ACC, TOF and CDC. 
Charged tracks from B meson candidates 
with TZk > 0.6 are identified as kaons and 
TZk < 0.4 as pions. Kaons produced in the 
B ± — ► ujK decays are selected by these cri- 
teria with efficiencies 85%, while the corre- 
sponding rates of Kaons that are misiden- 
tified as pions are 11%; on the other hand, 
pions from B^ — > ujtt^ are selected with an 
efficiency of 90% and a kaon misidentifica- 
tion rate of 9%. Candidate 7r° mesons are 
reconstructed from pairs of photons with in- 
variant mass in the range 0.1178 GeV/c 2 < 
M 77 < 0.1502 GeV/c 2 . Candidate K° s 
mesons are reconstructed using pairs of op- 
positely charged particles that have an in- 
variant mass in the range 0.482 GeV/c 2 < 
M n+n - < 0.514 GeV/c 2 . The vertex of the 
candidate is required to be well recon- 
structed and displaced from the interaction 
point, and the Kg momentum direction must 
be consistent with the K s flight direction. 
Candidate u — > 7r + 7r~7r° decays are recon- 
structed from charged tracks with TZk < 0.9, 
where the pions from the u decay are se- 
lected with 96% efficiency, and the 7r°s with 
center-of-mass frame momentum are greater 
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than 0.35 GeV/c 2 . Candidate uj mesons 
are required to have invariant masses within 
±30 MeV/c 2 of the nominal value. 

B meson candidates are formed by com- 
bining an uo meson with either a kaon 
(K ± , Kg) or a pion (it^, 7r°). Two kine- 
matic variables are used to select B candi- 
dates: the beam constrained mass M^ c = 
\/(£bcam) 2 - I^b M | 2 and the energy differ- 
ence AE = E™ - E™ m , where E™ m is 
the beam energy in the T(4S) rest frame, 
and Pb M , Eg M are the momentum and en- 
ergy of the B candidate in the T(4S) rest 
frame. Candidates with M^ c > 5.2 GeV/c 2 
and \AE\ < 0.25 GeV (\AE\ < 0.30 GeV for 
B° — > uj7t°) are selected for further analysis. 

The dominant background comes from 
quark- ant iquark continuum events (e + e~ — > 
qq, q = u,d,s,c). The continuum background 
is characterized by a jet-like structure while 
the BB events have a more spherical dis- 
tribution. Several event-shape variables are 
employed to suppress the continuum. The 
thrust angle 6t is defined as the angle be- 
tween the thrust axis || of the B candidate 
daughter particles and that of the rest of 
the particles in an event. Signal events are 
uniformly distributed in cos 8t, while contin- 
uum events are sharply peaked near cos 6? = 
±1.0. Events with |cos#t| < 0.9 are se- 
lected. A Fisher discriminant is formed by 
combining a set of modified Fox- Wolfram mo- 
ments with the variable S±, defined as the 
scalar sum of the transverse momentum of 
all particles outside a 45° cone around the 
thrust axis of the B candidate divided by the 
scalar sum of their momenta. Further vari- 
ables that have been found to separate signal 
from continuum background include: the co- 
sine of the angle between the flight direction 
of the B candidate and the beam direction 
(cos^b), the distance along the beam direc- 
tion between the B vertex and the vertex of 
the remaining particles of the event (Az) and 
the cosine of the helicity angle, defined as 
the angle between the B flight direction in 
the to rest frame and the normal direction 



of the plane spanned by the three daugh- 
ter pions of the uo. The probability density 
functions (PDFs) for these three variables 
and the Fisher discriminant are obtained us- 
ing events in signal Monte Carlo simulation 
and data with Mb c < 5.26 GeV/c 2 for signal 
and qq background, respectively. These vari- 
ables are combined to form a likelihood ratio 
1Z = C J^ S C where £s(qq) is the product of 
signal (qq) probability densities. 

Additional background discrimination is 
provided by the quality of the B flavor tag- 
ging of the accompanying B meson. We use 
the standard Belle B tagging package (loT |. 
which gives two outputs: a discrete variable 
(q) indicating the B flavor and a dilution fac- 
tor (r) ranging from zero for no flavor in- 
formation to unity for unambiguous flavor 
assignment. We divide the data into six r 
regions. Continuum suppression is achieved 
by applying a mode dependent requirement 
on TZ for events in each r region based on 
M s I \JM S + M q q } where Jsf a is the expected 
number of signal events estimated from sim- 
ulation and M q q denotes the number of back- 
ground events estimated from data. This 1Z 
requirement retains 74%, 68%, 74%, and 57% 
of the signal while rejecting 91%, 94%, 91%, 
and 95% of the continuum background for 
the uK ± , a;7r ± , uK°, and un° modes, respec- 
tively. 

Simulation studies indicate small back- 
grounds from generic b — > c transitions in 
the charged B modes; they are found to be 
negligible for the neutral modes. Moreover, 
two other backgrounds need to be consid- 
ered: the ujK + <r-* um + reflection, due to 
K + <r-* 7r + misidentification, and the feed- 
down from charmless B decays, predomi- 
nantly B -> u K* (892) and B -> up(770). 
We include these three components in the fit 
used to extract the signal. 

The signal yields and partial rate asym- 
metries are obtained using an extended un- 
binned maximum likelihood fit on M^ c and 
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AE. The likelihood is defined as 



(2) 



i j 



Vj = -[l-Qi-AcPjjPjiM^AEi) (3) 

where % is the identifier of the z-th event, 
P(M bc , AE) is the PDF of M bc and AE, Q 
indicates the B meson charge, B + (Q = +1) 
or B~{Q = — 1), Mj is the number of events 
for the category j, which corresponds to ei- 
ther signal, qq continuum, a reflection due to 
K-7T misidentification, or BB backgrounds, 
inclusive of b — > c background and other 
charmless -B decays. For the neutral B mode, 
Vj in Eq.(2) is simply Pj(M hci , AEi) and 
there is no reflection component and b — > c 
background. 

The signal distribution in Mb c is parame- 
terized by a Gaussian function centered near 
the mass of the B meson, while the Crystal- 
Ball line shape llj is used to model the AE 
distribution. The product of the two func- 
tions is used as the signal PDF, calibrated 
with large control samples of B — > Dtt, Z5° — > 
K + 7i^7T and B — > i]'(i]'n- + 7T~)K decays. The 
continuum PDF is the product of a first or- 
der polynomial for AE and an ARGUS func- 
tion |l2| for Mb c . Other background PDFs 
are modelled by a smoothed two-dimensional 
Mb c — AE function obtained from the MC 
simulation. The yields of the BB back- 
ground mentioned in the previous paragraph 
are fixed with Acp = in the fit based on 
the MC simulation except for the reflection 
component. The Acp and the normaliza- 
tions of the reflection components are fixed 
to expectations based on the B + — > uiK + 
and B + — > uoti + partial rate asymmetries and 
branching fractions, as well as K + <-> 7r + fake 
rates. The reflection yield and Acp are first 
input with the assumed values and are then 
recalculated according to our measured re- 
sults. 

Table |l] and Table [H] show the measured 
branching fractions and Acp as well as other 
quantities associated with the measurements. 



Branching fractions are computed as the sum 
of the yields divided by the corresponding ef- 
ficiencies in each dataset, divided by the to- 
tal number of B mesons. The reconstruction 
efficiency for each mode is defined as the frac- 
tion of the signal yield remaining after all se- 
lection criteria, where the yield is determined 
by performing the unbinned maximum like- 
lihood fit on the MC sample. We define the 
significance as J— 21n(£ /£ max ), where £ max 
is the maximum likelihood from the fit when 
the signal branching fraction is floated, and 
Co is the likelihood obtained when the signal 
branching fraction is set to zero. Systematic 
uncertainties are included by repeating the fit 
with all parameters fixed at ±lcr from their 
central value and choosing the set of param- 
eters that gives the smallest significance. 

For the branching fraction measurement, 
the main sources of systematic uncertainties 
are: charged tracking efficiency (1.0% per 
track); neutral pion detection (3.5% for Set I 
and 4.0% for Set II); K s reconstruction (4.5% 
for Set I and 4.0% for Set II); to mass reso- 
lution (0.2% for Set I and 0.6% for Set II); 
requirement on 1Z (about 3.0%). System- 
atic uncertainties due to fit parameters are 
evaluated by varying each of the fixed pa- 
rameters by ±lcr, and by adding in quadra- 
ture the resulting deviations from the central 
value. The total systematic uncertainty is the 
quadratic sum of all the above contributions. 

In summary, we update the measurements 
of the branching fractions of the exclusive 
two-body charmless hadronic B decays with 
an uo meson in the final state. We measure 
the following branching fractions: B{B ± — > 
uK*) = (8.1 ± 0.6 ± 0.5) x 10" 6 , B{B ± -> 
wtt*) = (7.0 ±0.6 ±0.5) x 10" 6 , and B(B° -> 
wK ) = (3.9 ±0.7 ±0.4) x 10- 6 . We also set 
the following upper limit: B(B° — > urn ) < 
1.5 x 10~ 6 at the 90% confidence level. Fi- 
nally, we measure the partial rate asym- 
metries in B^ — > ujK^ and 5 ± — > ojh 
A CP = 0.05 ± 0.08 ± 0.01 for uK*, and 
A CP = -0.03 ± 0.09 ± 0.02 for cjtt*; no ev- 
idences of direct CP violation is found in ei- 
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TABLE I: Signal yields (M s ), efficiencies (etot), 
fitted branching fraction (B), partial rate asym- 
metry (Acp), and fit significance including sys- 



tematic error (£) for the B~ 



UJ 



K± and ir 



B- 



ujit 



modes. 



Mode 






M s (Setl) 






M a (Setll) 


170 ± 16 


145 ± 16 


e to t (Setl) 


8.6 ±0.1 % 


9.1 ±0.1 % 


e to t (Setll) 


8.5 ±0.1 % 


8.6 ±0.1 % 


B (xHT 6 ) 


8.1 ±0.6 ±0.5 


7.0 ±0.6 ±0.5 




0.05 ±0.08 ±0.01 


-0.03 ±0.09 ±0.02 


£ 


20.0o- 


17.3(7 




AE (GeV) 




AE (GeV) 



(GeV/cO 



TABLE II: Signal yields (M s ), efficiencies (etot), 
fitted branching fraction (B), upper limit on the 
branching fraction (U.L.) including systematic 
error, and fit significance including systematic 
error (£) for the B° — > luK° and B° — > U7r° 



modes. 



Mode 






Ms (Setl) 


16+ 5 - 
10_ 44 




A4 (Setll) 


1 4 +5.1 


2.811:1 


e tot (Setl) 


2.7 ±0.1 % 


4.7 ±0.1 % 


e t ot (Setll) 


2.8 ±0.1 % 


3.7 ±0.1 % 


B (xl(T 6 ) 


3.9 ±0.7 ±0.4 


0.8±g;]| ± 0.1 


U.L. (xl0~ 6 ) 




1.5 


£ 


9.3cr 


2.2cj 



ther mode so far. 
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FIG. 1: Projections of fit results on Mt, c (within 
the region of \AE\ < 0.1, -0.15 < AE < 0.1 
for urn ) and AE (within the region of 5.27 < 
M hc < 5.29) for uK±, utt^, loK°, and uir°. 
Open histograms are data, solid curves are the 
fit functions, dash-dotted lines represent signals, 
the dashed lines show the qq continuum, light 
gray crossed-hatched areas represent rare B de- 
cay backgrounds, dark gray hatched areas peak- 
ing in the M^ c signal region represent reflections 
due to K-tv misidentification and the small dark- 
gray contributions in the negative AE region of 
the charged modes are due to generic BB back- 
ground. 
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FIG. 2: Projections of fit results on Mb c (within 
the region of \AE\ < 0.1, -0.15 < AE < 0.1 
for ton ) and AE (within the region of 5.27 < 
M bc < 5.29) for wK*, um*, loK°, and uir°. 
Open histograms are data, solid curves are the 
fit functions, dash-dotted lines represent signals, 
the dashed lines show the qq continuum, light 
gray crossed-hatched areas represent rare B de- 
cay backgrounds, dark gray hatched areas peak- 
ing in the Mb c signal region represent reflections 
due to K-ir misidentification and the small dark- 
gray contributions in the negative AE region are 
due to generic BB background. 
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